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Abstract 



We discuss the contributions of higher order terms in weak Hamiltonian to the standard two-neutrino 
double beta decay. The formalism for the unique first forbidden transitions has been developed, and 
it is shown that they can alter the two-electron energy spectrum. Yet, their effect is too small to 
screen the detection of exotic neutrinoless double beta decays, which are candidates for testing the 
physics beyond the standard model. 
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The double beta {(3(3) experiments furnish a unique window onto whatever new physics may 
replace the Standard Electroweak Model (SM). To be observable the new physics should: i) violate 
the electron-lepton-number (L e ) conservation, that is fulfilled in the SM, and/or ii) fit scalar particles 
(currently called Majorons), light enough to be produced in the (3(3 decay. 

The quantity that is used to discern experimentally between the ordinary SM two-neutrino decays 
(P/hu) and the exotic neutrinoless (3(3 events, both without {(3(3q v ) and with Majoron emissions ((3(3 M ), 
is the electron energy spectrum dT/de of the decay rate T, as a function of the sum e = t\ + e 2 of the 
energies of the two emitted electrons. The (3(3 2v decay exhibit a continuous spectra in the interval 
2 < e < Q, where Q = E I — E F is the released energy. On the other hand, when no new light particles 
are created, the L e -violating terms in the weak Lagrangian, that generate a Majorana mass for the 
neutrino, can be identified if they produce the (3(3q u decay, with the energy spectrum that is just a 
spark at the energy Q. The spectra for the (3(3 M decays are also continuous, but their shapes are 
clearly differentiated from that for the (3(3 2v decay, and depend on whether one or two Majorons are 
emitted and on the leptonic charge (L e = 0, —1, —2) they carry [I], ^ [|. Thus, both (3(3 Qv and (3(3 M 
processes, that are potentially capable to reveal the new physics, are clearly distinguishable from the 
SM (3(3 2v decay. 

The sensitivity of the (3(3 decay experiments is steadily and constantly increasing. For instance, 
while the pioneer laboratory measurement of the (3(3 2v decay in 82 Se has relied on only 40 events [Q, 
the most recent experiment on 76 Ge |J was done with high statistics (~ 20000 counts). Another 
example is the evolution of the half life limit for the (3(3q v decay in 7G Ge. ^From the first measurement 
in 1952, T 0i , > 2 • 10 16 y @, it has varied to T 0v > 2.2 • 10 22 y in 1983 0, while the most recent 
value is T 0u > 1.2 • 10 25 y ||. By comparing the last one, as well as the measured half life limit for 
(3(3 M decay: T M > 1.67 • 10 22 y, with the corresponding half life for f3(3 2v decay: T 2u = 1.77 • 10 21 y, it 
can be said that presently are being observed effects of the order of 10~ 4 at e ~ Q and of the order 
of 10 _1 at e ~ Q/2. There are also several ongoing and planned experiments that are supposed to 
allow for measurements of still smaller effects. The most promising one seems the GENIUS project 
f|, which is supposed to test the (3(3q v half-life of 76 Ge up to a limit of T 0u > 5.7 ■ 10 28 y. One might 
expect that the sensitivity for probing the (3(3 M decay will be improved accordingly as well. 

The measured limits on the exotic (3(3 transition probabilities are being rapidly translated into 
more stringent constrains on the parameters of new theoretical developments in particle physics, 
such as: Majorana mass of the light neutrinos, right-handed weak couplings, right-handed weak 
coupling involving heavy Majorana neutrinos, massless Majoron emission, R-parity breaking in the 
supersymmetric model, etc. They have also broad consequences on the history of the primordial 
universe, evolution of stellar object and astrophysics of supernovas. 

In confronting the experimental data with theory, the validity of allowed approximation (A) is 
usually assumed for the standard (3(3 2v decay. This implies to consider only virtual states with 
spin and parity J w = + and 1 + , which contribute via the nuclear operators g v ^ + and g A T + er, 
respectively. The higher order effects, coming from the parity-forbidden (PF) virtual states with 
J 71 " = 0~, 1~, 2~, have been ignoring almost entirely all along by workers in the field, simply because 
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they expected them to be small. But, in planing future experiments and in searching for exotic decays, 
it might be important to know how small these effects are, and whether they could eventually lead 
to experimental consequences similar to that produced by the former. 

Besides, in the charged Majoron (CM) model designed by Burgess and Cline [Q and by Carone 
0, which is the most hopeful one to be observed experimentally among the new Majoron models ||, 
the (3(3 decay proceeds via relativistic corrections in the hadronic current. More, the nuclear matrix 
element is of the form M. CM = M.~^ M — -M-cmi with A4^m De i n g the contributions of two heavy Dirac 
neutrinos with masses M±, and there is a strong destructive interference between M.^ M and M.~^ M 
when M + = M_ || . Therefore it might be also interesting to compare the outcome of the CM model 
with that arising from the PF transitions in the standard (3j3 decay. [] 

Naively thinking it could be inferred that the PF contributions to the (3(32 V amplitude are of 
the order of (RQ/4) 2 , where R is the nuclear radius and the factor 4 comes from the fact that the 
Q- value is shared by 4 leptons. Therefore, as for medium heavy nuclei (RQ/4) 2 ~ 10~ 3 , they would 
alter the half lives and spectrum shapes at the level of 10~ 3 or 10 -6 , depending on whether there are 
interference or not between the A and PF contributions, i.e. whether the PF matrix elements enter 
linearly or only quadratically in the expression for the decay rate. 

The above estimate, however, it is not valid for the non-unique (NU) transitions with AJ n = 
0~, 1~, due to: i) the electron s- waves contributions via the velocity dependent terms in the hadronic 
current, and ii) the Coulomb enhancement of the matrix elements coming from the pi/2 waves. When 
these effects are considered, and because the interference between the NU and A matrix elements, the 
theoretical half lives are increased in ~ 30% [[□], |12[. Yet, due to the nuclear structure uncertainties, 



we cannot disentangle from the measured half lives alone how large is the effect of the NU matrix 
elements, even when the experimental errors are relatively small. 

We have also pointed out |Ll| that the NU transitions do not modify the allowed shape of the 



two-electron spectrum, and that the unique (U) transitions states should be examined. It should be 
remembered that the spectrum shape of a single (3 transition of this type, provides for the emission of 
more high- and low-energy electrons, than are found in spectra that have the allowed shape |13|, |14| . 
Thus, it sounds appropriate to speculate about a similar effect in the (3(3 decay. Simultaneously, it is 
essential to find out whether the U matrix elements contribute linearly or quadratically to the (3(3 2l , 
decay rate. 

Meanwhile, Civitarese and Suhonen ||15|| , have pointed out that the effect of the U transitions on 
the T 2u in 76 Ge can be disregarded, but they have not discussed at all the corresponding spectrum 
shape. Without performing any theoretical development, it is simply assumed that the U matrix 
element come through quadratically, and that the ratio between the phase-space factors, for A and 
U (3(3 2v transitions, is of the order of 10 6 . Besides, in the same work it is suggested that the NU 
transitions should be retarded as well by the factor (RQ/4) 4 ~ 1CU 6 . 

Before presenting the numerical results for the [3 feu decay through the virtual states J 71 = 2~ , 

1 In fact, except for the neutrino mass term, the (5(3® v decay also arises from the higher order effects [jlO|. 
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we briefly sketch the derivation of the corresponding formulae, which has not been done so far. The 
/302V decay rate reads 



dY 2v = 2n^\R 2l/ \ 2 5(ti + e 2 + 0J1 + uj 2 - Q) f[ dp k dq k 



k=l 



where the symbol represents both the summation on lepton spins, and the integration on neutrino 
momenta and electron directions. For a transition from the initial state |0 7 ) in the (N,Z) nucleus 
to the final state \0 F ) in the (JV — 2, Z + 2) nucleus (with energies E l and E F and spins and parities 
J 71 = + ) the transition amplitude is evaluated via the second order Fermi golden rule: 

A ZlT ) n -fcjv — -C// + ei + 

where e« = (e,, pj, s ej ), U{ = (u{, q^, s Ui ), P{l\l 2 ) exchanges the quantum numbers of leptons l\ and l 2 , 
and n runs over all levels in the (JV — 1, Z + 1) nucleus. The weak Hamiltonian reads 

# w (ei/) = ^= / dx^(x) J^(x) + /i.e., (3) 

where G = (2.996±0.002)xl0 -12 is the Fermi coupling constant (in natural units), j M (x) is the usual 
left-handed leptonic current , and for the hadronic current 

J"(x) = (pv(x)-p A (x),j v (x)-j A (x)), (4) 

the following non-relativistic approximation will be used [ffO| 

/EV(x) = £>^r+5(x- r n ), 



Pa(x) = T n [<*n ■ Pn^(x - T n ) + <J(x - r n )<T n ■ Pr 

JV, 

j A (x) = 5f A ^r+cr n 5(x- r n ), (5) 



2M N n 

jv(x) = ^^H r « [Pn5(x-r n ) + 5(x - r n )p„ + / w Vx(T„J(x - r n )] 

Z1V1 N n 



where M N is nucleon mass, and g v , g A and f w are, respectively, the vector, axial- vector and weak- 
magnetism effective coupling constants. 

In the discussion of the f3(3 2v decay we ignore both the weak-magnetism term, and the action 
of the velocity dependent terms on the lepton current. These terms cause the "second-forbidden" 
contributions, which do not alter the electron spectrum shape and will be discussed elsewhere. Ad- 
ditionally, it will be assumed that the Coulomb energy of the electron at the nuclear radius is larger 
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than its total energy, which leads to the ^-approximation fll2| , |14f . Thus, for the purposes of the 
present study, and after a lengthy algebra, we cast the weak Hamiltonian in a rather novel form 



G 



(6) 



which can be used for multiple purposes. Here Wj and Wj are, respectively, the allowed and 
forbidden nuclear operators, and 



LAeu) 



sg{s l 



2e 



v ■ P 

e + 



(7) 



are the leptonic matrix elements, with xi s ) being the usual Pauli spinor. The leptonic operators lj 
are listed in Table [l], together with WJ. 



Table 1: Operators tj and WJ for different multipoles J; p = p[Fi(e)/F (e)] 1//2 , v = p/M N and 
£ = aZ/2R. 



J 


ij 


w+ 


Wj 







1 


9v 


-9a(<t ■ v + 


£i<x • r) 


1 


a 




-fiW - £[g v ir ■ 


-g A (crxr)} 


2 [<r<g 


Kq + p)] 2 




ig A (a ® r 


) 2 /\/5 



In the next step we evaluate the transition amplitude and get 
G 2 



4(2tt) ( 



L (en/i) • L (e 2 ^ 2 ) + M%j) - L 2 ( ei ^) • L 2 (e 2 z/ 2 )A4 



where -M 2l , 



AM0+) + M 2i ,(l+), M£f = M 2 ,(0-)+M 2 ,(l-), and Mg, = M 2v {2~) are, respec- 
tively, the (3f32v matrix elements for the A, NU and U transitions, and 



Of) 



(9) 



After introducing (^) into ([!]) and performing the spin summations and angular integration, 
the contribution of the lepton matrix elements L 2 (ei^i) • l^^z^L^eii/j) • l_Q(e 2 z/j) turns out to be 
identically null for i,j = 1,2 or 2, 1. Therefore, there is no interference term between the A and U 
matrix elements., as happens with M. 2 v an d We get 



dT^ = dT^ NU + dT u 



2// 



4G 4 
15^ 



\m& + mi u \ 2 dni + \M u 2 j 2 dn u 2u 



(10) 
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where 

1 2 
dttl = 7^(Q-^- e 2) 5 II PktkF Q (e k )de k , (11) 
1 n k=i 

is the usual phase space for the f3(3 2 u in the A approximation, and 

dn 2» = W^TT E^(Q-ei - e 2 ) 5+2j f[ p k e k F (e k )de kl (12) 

with ao = P1P2, a i = 16p 2 /35 and a 2 = 1/21. 
The corresponding half life is 

^(OMC) = (S&\M& + Mg\ 2 + gi\MZ,\ a y\ (13) 

where 

4G 4 /■ 

2v 15tt 5 In 2 J 2u ' V ; 

are the kinematical factors. 

The spectrum shapes dT^ 1 * 11 /de and dT^/de for 76 Ge are confronted in Fig. 1. At variance with 
the single f3 emission, the spectrum for the U double beta process deviates from the allowed shape in 
the low-energy region, but not for e = Q. Thus, independently of magnitude of T 2u , the virtual states 
J 77 = 2~ will never interfere with the detection of the /3/3q u events. Their spectra, still, can overlap 
with those engendered by the (3(3 M decays. This is illustrated in same figure for the case of the CM 
model. 



Table 2: Kinematical factors Q 2u , and the nuclear matrix elements M. 2v evaluated within the QRPA 
formalism. 



Nucleus Gt \v~ x \ SI [y 



,-1 



M A 



76 Ge 5.39 10" 20 2.10 10" 19 
82 Se 1.80 10~ 18 2.54 10~ 17 
100 Mo 3.91 10- 18 5.50 10~ 17 



2u 



l 2u 



2v 



0.050 
0.060 
0.051 



-0.008 1.0 10" 5 
-0.009 9.8 10~ 6 
-0.014 1.1 10~ 5 



Numerical results for the kinematical factors and the nuclear matrix elements, for several exper- 



imentally interested nuclei, are displayed in Table 

pn-QRPA model, following the procedure adopted in our previous works 
seen that: 



The moments M. 2v were evaluated within the 

12 , 17] . It can be easily 



\M 



U I 
2v\ 



R 



+ Ml u \ 



$2u 



(Q/4) 



ir>A 
V2w 
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Therefore, from the theoretical developments and numerical calculations done here, it can be stated 
that the simple estimate 

K/K +NU = (RQ/4)- 4 ~ 10 6 , 
is appropriate for the unique transitions. 

Table 3: Calculated half-lives (in units of y) for the A+NU and U transitions and for the charge 
Majoron emission. 



Nucleus 


rpA+F 

l 2v 


rpU 

1 2u 




T cm (m+^M-) 


w Ge 


1.1 1(F 


4.7 10 28 


1.6 10 25 


1.3 10 2y 


82 Se 


2.1 10 20 


4.1 10 26 


9.2 10 23 


9.5 10 27 


100 Mo 


1.9 10 20 


1.4 10 26 


3.8 10 23 


4.1 10 27 



Finally, in Table [| are compared the half-lives for the /3/3z v and f3f3 CM decays. In the CM model 
the effective coupling constant was taken to be g CM = 9 2 /2 with 9 = 0.1, and two different values 
for the heavy Dirac neutrinos masses were considered, namely M + — > oo and M = 100 MeV, and 
M± = MVT±9, with M = 100 MeV |, ||. It turns out that: 

T cm {m + ->oc)/T£ nu ~ 10 3 ; T cm (m + ^)/T 2 1 +nu ~ 10 7 . 

Thus, the emission rate for the recently discovered Majoron models is very strongly conditioned by 
the model parameters, and it could be so small as that arising from the unique forbidden transitions. 

In summary, the higher order effects in the standard physics modify the flfyv spectrum shape 
but only at the level of 10~ 6 and mainly at low two-electron energy, where most backgrounds tend 
to dominate. Therefore they would hardly mask the observation of the potential exotic (3f3 decays. 
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Figure Captions 

Fig. 1. Electron energy spectrum for the nucleus 76 Ge, as a function of the sum of energies of the 
two emitted electrons, for: the standard 2v allowed (PP^u) an d unique forbidden transitions 
and the exotic neutrinoless decays, with Majoron charged emission {(3(5 CM ) and without (/3j3ou)- All 
four curves have been arbitrarily assigned the same maximal values for purposes of comparison. 
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